Adiabaticity and diabaticity in strong-field ionization 
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If the photon energy is much less than the electron binding energy, ionization of an atom by 
a strong optical field is often described in terms of electron tunneling through the potential bar- 
rier resulting from the superposition of the atomic potential and the potential associated with the 
instantaneous electric component of the optical field. In the strict tunneling regime, the electron 
response to the optical field is said to be adiabatic, and nonadiabatic effects are assumed to be 
negligible. Here, we investigate to what degree this terminology is consistent with a language based 
on the so-called adiabatic representation. This representation is commonly used in various fields of 
physics. For electronically bound states, the adiabatic representation yields discrete potential energy 
curves that are connected by nonadiabatic transitions. When applying the adiabatic representation 
to optical strong-field ionization, a conceptual challenge is that the eigenstates of the instantaneous 
Hamiltonian form a continuum; i.e., there are no discrete adiabatic states. This difficulty can be 
overcome by applying an analytic-continuation technique. In this way, we obtain a rigorous classifi- 
cation of adiabatic states and a clear characterization of (non)adiabatic and (non)diabatic ionization 
dynamics. Moreover, we distinguish two different regimes within tunneling ionization and explain 
the dependence of the ionization probability on the pulse envelope. 

PACS numbers: 32.80.Rm, 42.50.Hz, 31.15.A- 
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I. INTRODUCTION 



The realm of strong-field physics has become a focal 
point of interest in the atomic, molecular, and optical 
physics community over the last two decades. This was 
particularly supported by the rapid development of lasers 
producing high intensities (10 14 — 10 15 W/cm 2 ) that gen- 
erate forces comparable to intra-atomic forces, and ul- 
trashort pulse durations of the order of femtoseconds 
(lCT 15 s) down to attoseconds (lfr 18 s) [Hi]. The time- 
resolved investigation of electron dynamics in atoms and 
molecules has come into reach because the typical time 
scales involved in electronic excitations (between 50 as 
and 50 fs) can be accessed. 

The process of tunneling ionization has been studied 
extensively. Following the calculation of the tunneling 
ionization rate for the ground state of hydrogen in a 
static electric field by Landau Q, Keldysh extended the 
theory to ionization by strong electromagnetic fields |J]. 
Later, Ammosov, Delone, and Krainov ("ADK") general- 
ized the results to slowly varying fields by introducing the 
quasistatic approximation and defining the tunneling ion- 
ization rate by averaging over one optical period ( "ADK 
theory" ) p . A self-contained derivation of the tunneling 
rate in this approximation is presented in Ref. |6|. In the 
original derivation |^| Keldysh introduced the parame- 
ter 7 = y/IrJ(2U p ), which is now known as the Keldysh 
parameter [7|. Here, I p is the ionization potential and 
Up is the ponderomotive potential, which corresponds to 
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the average energy of a free electron oscillating in the 
electric field. According to Keldysh, 7 divides the phe- 
nomenon of strong-field ionization into two regimes: for 
7 <C 1 ionization is governed by tunneling ionization [3[ , 
while for 7>1 the process is governed by perturbative 
multiphoton ionization [8J. In the range of 7 « 1 both 
effects compete with each other [9|, [HJ ■ In later papers 
the Keldysh parameter has been connected to the notion 
of adiabaticity of the ionization process [n], [l~2| ■ Far into 
the tunneling regime, the atomic response is considered 
to be purely adiabatic. Adiabatic means in this context 
that the ionization rate at a given time is solely defined 
by the instantaneous electric field. 

More generally, when a time-dependent process is adia- 
batic, the state of the system at any given time is always 
an eigenstate of the instantaneous Hamiltonian, which 
depends on one or more external parameters (like the 
electric field). Consequently, the energy eigenstates and 
their corresponding eigenenergies become parametrized 
and lead to energy curves (or energy hypcrplanes de- 
pending on the number of external parameters). Nona- 
diabatic dynamics occur when transitions between adi- 
abatic curves start to appear. This is, particularly, the 
case when two adiabatic curves are energetically close to 
each other and the external parameters are changed rela- 
tively fast such that the system has no time to "instanta- 
neously" respond to the change. As a result, the system is 
not in one defined adiabatic state anymore but rather in a 
superposition of several adiabatic eigenstates. In various 
fields of physics and chemistry the adiabatic representa- 
tion has been used to study adiabatic and nonadiabatic 
effects. Its application includes fields like Rydberg atoms 
J13l Il4j , molecular dynamics |15l - fl7j , atomic and molec- 
ular collisions |18l - (20l |. and ultracold gases and trapped 



ions [21H23|. 

An important aspect in the adiabatic representation is 
the discreteness of eigenstates which is essential to obtain 
a discrete set of energy curves. In the case of strong-field 
ionization, however, the instantaneous eigenstates of the 
Hamiltonian form a continuum. Therefore, the identifi- 
cation of a nonadiabatic effect happens rather indirectly 
[24J |25| : either the spectrum of the photoelectron after 
the pulse or the field dependence of the ionization rate 
is analyzed. Various results on nonadiabatic behavior 
in stron g-fie ld ionization have been presented in the lit- 
erature [26H28J and there are many different usages of 
the terms "adiabatic" and "nonadiabatic" . By introduc- 
ing an analytic continuation in the complex plane the 
instantaneous Hamiltonian becomes non-Hcrmitian and 
tunneling states appear as discrete eigenstates. These 
discrete states can be now used to apply the adiabatic 
representation to strong-field ionization dynamics. 

In this paper we strictly apply the adiabatic repre- 
sentation to strong-field ionization and find that in the 
tunneling regime the ionization dynamics is defined by a 
diabatic rather than an adiabatic behavior. Diabatic dy- 
namics means that the response of the system follows one 
specific diabatic state. Here, the diabatic states are de- 
fined by the overlap with the field-free eigenstates. In this 
formulation we find that the ionization dynamics can be 
divided into two regimes. Furthermore, with increasing 
frequency we observe a transition from the diabatic to the 
nondiabatic regime. In particular, we study the few-cycle 
limit and find a non-constant population as a function of 
the optical frequency which has been interpreted in the 
literature as a sign of a nonadiabatic process [1J, |25| . 
We show for a few-cycle pulse with a Kcldysh parameter 
7 < 1 that this effect rather represents a dependence 
on the form of the pulse and can be fully explained by 
a diabatic picture depending on a single diabatic state 
connected to the field-free ground state. The main text 
is divided into three sections: 

• Section [IT] is devoted to the general theory of the 
equations of motion in the adiabatic basis and in- 
troduces also diabatic states. 

• The third section presents one-photon absorption 
as an extreme case of a nonadiabatic/nondiabatic 
ionization process. 



• In Sec. IIV1 the central section of this paper, we 
develop the concept of diabaticity in strong-field 
ionization. We examine the transition from the di- 
abatic to the nondiabatic ionization regime. 

Atomic units are employed throughout unless otherwise 
indicated. 



II. ADIABATIC EIGENSTATES 



adiabatic. In the following, we derive the quantum- 
mechanical equations of motion in the adiabatic basis, 
which is given by the states that are eigensolutions to 
the Hamiltonian of the system for a set of instantaneous 
parameters. 

Let us study a system where the Hamiltonian depends 
on an external time-dependent parameter e(t). The time- 
dependent Schrodinger equation has the form 

itkmt)) = H(t)Mt)) = {H, + U[e(t)}} |*(i)>. (1) 

Hq describes the atomic Hamiltonian, whereas U includes 
all external potentials and is dependent on the parameter 
e(t). At a given time t, the instantaneous eigenstates, 
which constitute the adiabatic basis, are defined by 1 



H + U(t) (*„(£)> =£W(£)|tf„(£)). 



(2) 



To analyze adiabatic and nonadiabatic effects we expand 
the electronic wavefunction in terms of the adiabatic 
eigenstates, *(£)) = Sn^Wl^nW)- Upon inserting 
this expression into Eq. (TTJ and projecting onto the eigen- 
state |\I/ m (£)), the equation of motion for the coefficient 
a m (t) reads 



>(*)■ 



-iJ2 a n(t)(^r 



Xt)\d t \* n (t)) = a m (t)E m {t). (3) 



The off-diagonal matrix elements {9 m (t)\dt\^f n (i)) intro- 
duce couplings between different adiabatic eigenstates, 
thus making the dynamics nonadiabatic [29(. In the adi- 
abatic approximation, where these couplings are consid- 
ered to be very small, Eq. (J3J) becomes 

ia m (t) + »Om(*)<*m(*)|*m(*)> = Ot m (t)E m {t), (4) 

which is solved with the initial condition a m (0) — 1 by 



i(t) = exp 



-i / dt'E m (t') 



exp[«7 m (£)], (5) 



where 7 m (t) = i f* dt'($ m (t')\^m(t')), so that the sys- 
tem evolves in a specific adiabatic eigenstate with a 
phase. If, on the other hand, ( v I'm(£)|^'n(£)} cannot be 
neglected, the whole sum in Eq. <j3j> has to be considered, 
so that different adiabatic eigenstates get coupled and 
nonadiabatic motion emerges. We can use dt = |jf <9 e and 
express the off-diagonal coupling elements also in terms 
of the change in e: 



<*m|*n> = (*m|0 e *r 



de 



(6) 



Whenever a system is given time to adjust to the pa- 
rameters on which it depends, the response is called 



1 The time dependence is implicit via the parameter e(i). 



Considering a two-level system with an external per- 
turbation proportional to e, the Hamiltonian of the sys- 
tem takes the form 
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(7) 
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where A is an internal coupling parameter. In Fig. [T] the 
energy curves of the two adiabatic states |^i} and \^f 2 ) 
of this system are shown as a function of the external 
parameter e, assuming A = 1. The internal coupling be- 
tween the diabatic states |1) = (1,0) T and |2) = (0, 1) T 
results in the effect that the two adiabatic curves do not 
cross. This phenomenon is known as an "avoided cross- 
ing". We see that A is the energy splitting between |^i) 
and |^2) at the degeneracy point of the states |1) and 
1 2). If the parameter e is changed sufficiently slowly, the 
system will remain in a given adiabatic state |^j) if, for 
e < 1 or e > 1, the system was in the state V^i). Note 
that in the vicinity of e = 1 the character of the adiabatic 
states changes from |1) to |2) and vice versa. If e changes 
rapidly in the vicinity of e = 1, the system has no time 
to change the character of its state; it makes a transi- 
tion from one adiabatic state to the other and follows 
the diabatic states |1) and |2), respectively. These jumps 
between adiabatic curves make the resulting dynamics 
nonadiabatic. For a given value of the external parame- 
ter, we can obtain the diabatic states also by choosing the 
adiabatic eigenstates with the maximal overlap with the 
free states (e = 0). For e < 1, the diabatic state |1) has 
the maximal overlap with the adiabatic state |^i), while 
for e > 1 the overlap of state |1) with the adiabatic state 
^2) is maximal, and vice versa for the diabatic state 
1 2). Asymptotically, the states |1) and |2) correspond to 
the states |^i) and "J 2} before the avoided crossing, and 
vice versa after the crossing. Near the crossing an inter- 
polation is performed in order to obtain a continuous and 
smooth state. 

A system's dynamics can of course also be formulated 
in other representations, e.g., in a diabatic basis [30|, 
where the diabatic states do cross (see the states |1) and 
1 2) in Fig. [Tj). Usually the basis is chosen such, that the 
off-diagonal couplings in Eq. ([6]) vanish or are at least 
small Qjl, [3l|. However, the diabatic basis, which is 
derived from the adiabatic basis by a unitary transfor- 
mation, is not unique and there are many different ap- 
proaches for reaching a diabatic representation |32H34f . 
One practicable method of diabatization is a local dia- 
batization method, which means that the diabatic state 
is constructed pieccwisc in a two-level model: At each 
avoided crossing between two adiabatic states the dia- 
batic state is followed. To this end, the size of the over- 
lap with the corresponding field-free state can be used 
as a criterion. This method turns out to be fruitful for 
the description of diabatic and nondiabatic strong-field 
ionization (see Sec. II V[) . Once a diabatic representation 
has been found, one can ask with which rate transitions 
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FIG. 1: The energy curves of the two adiabatic states \^i) 
and 1^2} are shown as functions of the parameter e. Through 
the off-diagonal matrix element A/2 a nonadiabatic transition 
is possible, whereupon the system follows the diabatic states 
|1) and 1 2), respectively. 



between diabatic states occur. These transitions will be 
called nondiabatic. 

In the following section we will make use of the fact 
that for weak perturbations the adiabatic eigenstates 
can be approximated through the field-free eigenstates. 
Therefore, the diabatic states exhibiting the maximal 
overlap with the field-free states are also the adiabatic 
states. In this case, the nondiabatic transitions are ex- 
actly the nonadiabatic transitions described above. 



III. ONE-PHOTON ABSORPTION 

First, we analyze the case of one-photon absorption 
within the adiabatic representation. If the system is 
exposed to a weak electric field of the form F(t) = 
F cos(uit) (in the dipolc approximation, see Sec. HV[) . 
with a frequency w, the system Hamiltonian is perturbed 
by the term F(t) z [351 ] . where z is pointing in the direc- 
tion of the field (which is assumed to be linearly polar- 
ized). The Hamiltonian in Eq. ([I]) takes the form 



H(t) = H + F(t)z, 



(8) 



where Hq is the atomic Hamiltonian and the electric field 
F(t) is coupled classically to the dipole operator z of the 
electron [the field F(t) corresponds to the parameter e of 

Sec. in]. 

In the following, we show that in the adiabatic repre- 
sentation the off-diagonal coupling elements in Eq. ([3]) 
are crucial for introducing transitions. Let {^h, }^Lo be 
the eigenstates of the field- free Hamiltonian, Ho\^n ') = 

Wnl^n }. For simplicity, we assume that the initial and 
final states of interest in the one-photon transition are 
nondegenerate. Performing static perturbation theory to 



first order, the adiabatic eigenstates read [36 



In / 



, f) + E i!^l l <. t 



fc#n 



Wn — ^k 



(9) 



Inserting Eq. ([H]) in Eq. ^ with e being the field F, we 
obtain the nonadiabatic coupling elements to first order 
inF: 
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(10) 

We are now ready to solve Eq. (|3|) including nonadiabatic 
coupling. We may treat the operator Vp = C 9j? as a 
perturbing time-dependent operator and, hence, analyze 
the states with time-dependent perturbation theory j36| . 
The first-order correction to the zeroth-order coefficient 
[Eq. (J5j)] is given by 

Assuming ujj > u>i , we obtain the total transition proba- 
bility per unit time 



„W|2 






2-E 



:* 



(0) 



F z 
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vp 



(0), 



5(iOf— UJi— lo). 

(12) 



This equation is exactly Fermi's golden rule [37] . In the 
present approach it is the nonadiabatic coupling that in- 
duces one-photon transitions between the field-free eigen- 
states. Viewed in this way, the phenomenon of one- 
photon absorption is entirely nonadiabatic. In the one- 
photon case, the states are well separated by a large en- 
ergy gap and there is no avoided crossing due to the weak 
field, which is only a perturbation to the field-free states. 
Note that the adiabatic states coincide with the diabatic 
states in the weak-field limit. 



i.e., F(t)z. Furthermore, the form of the Hamiltonian 
[cf. Eq. ([8])] is a result of the dipolc approximation, which 
holds in our case, because the size of the system of inter- 
est (a few A) is much smaller than the wavelength of the 
light pulse (« 1 nm) gSEJ. 




FIG. 2: The pure Coulomb potential of the helium atom 
(solid red line) is tilted in the presence of the electric field 
(dashed green line). The dotted black line denotes the field- 
free ground-state energy. 

In order to describe strong-field ionization dynamics, 
the Schrddingcr equation of the atom exposed to the field 
has to be solved nonperturbatively because perturbation 
theory fails for these high field strengths. As shown in 
Sec.[TT]in the adiabatic case the system will follow a given 
adiabatic state without making any transition. However, 
in the presence of a static electric field, electronically 
bound states become tunneling states, which means that 
there is ionization via tunneling. 

In the following, we will study helium as a concrete ex- 
ample to illustrate tunneling ionization within the frame- 
work of the adiabatic representation. 



IV. STRONG-FIELD IONIZATION OF ATOMS 

While in the case of one-photon ionization the photon 
energy necessarily exceeds the ionization potential, we 
will now examine the situation where the atomic system 
is irradiated by an intense electric field F(i) with a low 
photon energy, i.e., many photons arc needed to ionize 
the atom. 

When applying a strong external field [sec Eq. (|5J)] the 
effective potential seen by an electron gets tilted (see 
Fig. [2]). Therefore, a barrier of finite height is created 
through which the electron can tunnel. (If the electric 
field is so strong that the electron's energy lies above 
the barrier, the electron can just leave the atom without 
tunneling. This effect is called above-barrier ionization 
[38l |39|.) This tunneling picture of a tilted potential re- 
lies on the length form of the light-matter interaction, 



Constructing adiabatic and diabatic states for 
helium 



As already discussed (see Sec.[l|, in strong-field ioniza- 
tion the spectrum forms a continuum where a direct ap- 
plication of the adiabatic representation is inconvenient. 
To overcome this problem, a rigorous analytical contin- 
uation of the Hamiltonian can be performed by rotating 
the electron coordinates about an angle into the com- 
plex plane; this procedure is called complex scaling |42| . 
Another way to generate discrete eigenstates is to add a 
complex absorbing potential (CAP) to the Hamiltonian 
[43| . It can be shown that the latter method, which is 
conceptually easier, is closely connected to the complex 
scaling approach [44| . The key idea here is that for every 
tunneling state, i.e., every adiabatic atomic state that 
allows the electron to tunnel through the field-induced 



barrier, there exists a discrete eigenstate — a so-called 
Gamow vector [45( or Siegert state [46| — of the instan- 
taneous Hamiltonian. A Siegert state is associated with 
a complex energy and lies outside the Hermitian domain 
of the Hamiltonian. In fact, the associated wavefunction 
is exponentially divergent for large distances from the 
atom. Complex scaling or the use of a CAP eliminates 
the divergent behavior and renders the tunneling wave- 
function square integrablc. Thus, by making the Hamil- 
tonian non-Hcrmitian, it becomes possible to calculate, 
within Hilbert space, the complex Siegert energies of tun- 
neling states. The imaginary part of the Siegert energy 
E provides the tunneling rate T of each Siegert state by 
the relation T = -2 Im(E) 0,0. 

In order to obtain the instantaneous eigenstates we 
solve Eq. ([2]) with the Hamiltonian in Eq. ([8]) including 
a CAP. This yields the adiabatic eigenstates and corre- 
sponding eigenenergies of the atom shown in Fig. [3l A 
more detailed description of the methods used is given 
in the Appendix [X] We observe many avoided crossings 
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FIG. 3: The real part of the energy of the first adiabatic 
eigenstates as a function of a static electric field. The inset 
magnifies avoided crossings for small electric fields. 

among the higher adiabatic eigenstates for field strengths 
in the range below 0.01 a.u. (1 a.u.= 5.14 x 10 9 V/cm), 
while the ground state energy does not change signifi- 
cantly. One might wonder whether for sufficiently slow 
ramping of the electric field the atom follows the adi- 
abatic ground state. Indeed, for field strengths up to 
0.02 a.u. the adiabatic ground-state energy seems to re- 
main constant. But we know that the electric field can 
mix a whole manifold of excited states into the field- 
free states. When this happens, the adiabatic ground 
state loses the character of the field-free ground state 
(cf. Fig. [TJ). Analyzing the avoided crossings involving 
the adiabatic ground state around the field strength of 
0.02 a.u., we find that the ramping of the field has to be 
so slow that it lies in the radio frequency regime. There- 
fore, the system does not follow the adiabatic ground 
state for the frequency range of light usually employed in 



experiments (typically around 800 nm, corresponding to 
4 x 10 14 Hz). 

The electronic state follows the instantaneous eigen- 
state that has the maximal overlap with the field-free 
ground state. This is exactly the diabatic behavior de- 
scribed in Sec. [TH where the electronic state jumps from 
one adiabatic state to the other, keeping its field-free 
character. Here, we employ the diabatization method 
already alluded to in Sec. [ill where we construct the di- 
abatic state 1 4" i (£)) from the adiabatic basis {| , Fn(i))} 
using the criterion of maximal overlap with the field-free 
state |*- 0) ), i.e., 



■Ad) 



*i (*)> = l*n(*)>, where 



(13) 



,(()) 



K*n(t)|*r>l>l(*m(*)l*r)l. Vm^n 



(Oh 



This can be done as long as there is one distinct adiabatic 
state with a prominent character of the corresponding 
field-free state, so that the (orthogonal) complement of 
adiabatic states which are mixed in is small and can be 
ignored. The procedure works in principle also for ex- 
cited states. However, for excited states the condition 
of a small admixture breaks down already at low field 
strengths, such that this construction method works best 
for the field-free ground state. The overlap of the corre- 
sponding diabatic state |^q ) with the field- free ground 
state |\&g ) is always larger than 90% for field strengths 
considered here (see Fig. [3b). Figures[3Ja) anddfb) show 

the real part of the energy and the tunneling rate of |*q ) 
as a function of the electric field. The shift of the real 
part of the energy is well approximated by a quadratic 
behavior; for low field strengths below 0.1 a.u. the pref- 
actor is in accordance with the literature value of the 
polarizability of the helium ground state [49l . l50j . As 
expected, the tunneling rate increases considerably for 
sufficiently high field strengths. For field strengths larger 
than 0.07 a.u. the ionization rate is well captured by the 
analytic expression derived in the tunneling limit of the 
strong- field approximation [7j. 

Studying the adiabatic eigenstates and the avoided 
crossings reveals the suitability of the diabatic state con- 
structed as shown above for the description of strong- field 
ionization. The advantage of the diabatic basis is that 
the system follows one single diabatic state, which gives 
a clear and intuitive picture for the explanation of the 
physics in the tunneling regime. 



B. Ionization dynamics 

So far, the analysis was performed for the spectrum of 
adiabatic eigenstates, i.e., for static electric fields. Now 
we introduce dynamics by considering a Gaussian pulse 
of the form 

F{t) = f(t) cos(wi) = F e-* 2 / 2r2 cos(wi), (14) 
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FIG. 4: (a) The real part of the energy of the diabatic state 
l^lrg ), and (b) its tunneling rate, F = — 2Im(_E), are shown 
as a function of the electric field. fr ^ ^^^ """■'<■" " f l,T/ '\ 
with the field-free ground state. 



(c) The overlap of |# aj ) 



where Fo is the peak strength of the electric field, r is 
connected to the full width of the pulse at half maximum 
by r 2 = FWHM 2 /(81n2), and ui is the field frequency. 

We want to calculate the ionization probability out 
of the diabatic state \^f ) when applying this pulse. 
Let us assume that we have found a diabatic basis in 
which this particular diabatic state can be described by 
a coefficient a . Then the exact wavefunction reads 

*(*) = T,i<*i £) (t) i &i 1) (t). In analogy to the case of the 
adiabatic representation, equations of motion can be ob- 
tained for the coefficients in the diabatic basis where now 
coupling elements between the diabatic states imply non- 
diabatic transitions [cf. Eq. fl^J)]. If, in a "diabatic ap- 
proximation", the nondiabatic transitions are neglected 
we obtain the following equation of motion for the coef- 
ficients: 



«?\t) 



E 



id) 



,(d) 



WW, 



(15) 



,(./) 



where T { is the ionization rate of the diabatic state i. 
From the ionization rate of our distinguished diabatic 



(d), 



S d ) I 



state its population evolution P '(t) = \a Q (t)\ during 
the pulse can be inferred. To this end, the equation of 
motion for the probability of remaining in this particular 
diabatic state is calculated (we omit indices for the sake 
of readability): 



dP _ d 
~dt ~ dt 



\a(t)\ 2 = a*(t)a(t) + a*(t)a(t). 



(16) 



which can be analytically solved by separation of vari- 
ables: 



P(i)=cxp|-| f dt'T[F(t')}^ 



(18) 



with the initial condition P(t = — oo) = 1. Note that the 
rate depends on the external field. Inserting the tunnel- 
ing rate of the diabatic state in Eq. (fT5|) we calculate the 
diabatic ionization dynamics. Thereby we observe how 
much is ionized out of \^ ). Deviations from Eq. (fT8|) 
in the population dynamics can be attributed to nondia- 
batic behavior, i.e., transitions to other diabatic states. 

The results for four selected photon energies are shown 
in Fig. [5] for an electric field amplitude of Fq = 0.25 a.u. 
The pulse duration is kept constant so that we can study 
the ionization regime from few- to multi-cycle pulses. 
The exact result refers to the numerical solution of the 
Schrodinger equation [see Eq. (JTj)], where all dynamics 
are included, while the calculation of the diabatic curve 
via Eq. (fT8|) involves only the diabatic state \^> '). The 
gray-shaded areas in the background indicate the pulse 
intensity. In the frequency range shown, the evolution of 
the ground state population is well described by consid- 
ering only the single diabatic state. For u> = 0.3 — 0.8 eV 
[see Fig. [5^i)-c)] the difference between the numerically 
exact and the diabatic calculation is insignificant, while 
for u — 1.5 cV [see Fig.[5ji)] the discrepancy between the 
two methods becomes more noticeable. This is exactly 
the difference which gives us a measure of nondiabatic- 
ity. To clarify this further, a comparison between the 
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FIG. 5: Comparison of the ground-state populations cal- 
culated via numerical solution of the Schrodinger equation 
and via the rate equation (|18l) for the distinguished diabatic 
state for four different photon energies. The pulse intensi- 
ties are highlighted in the background: the pulse amplitude is 
-Fo = 0.25 a.u., and the pulse duration is 400 a.u. (~ 10 fs). 



Inserting Eq. (fTBj) in this equation the following rate 
equation for the population is obtained (cf. Ref. [401 ) : 



P(t) = -T[F(t)} P(t), 



(17) 



two methods is shown in Fig. [6] for a peak field strength 
of 0.2 a.u. by depicting the populations [Fig. [5k.)] and 
the relative difference [Fig. [Sb)] between them after the 
end of the pulse. One can clearly see that for sufficiently 
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FIG. 6: (a) Ground-state population after the end of the pulse 
calculated via numerical solution of the Schrodinger equation 
and from the single diabatic ground state as a function of 
the photon energy, and (b) relative difference between the 
two results, corresponding to the degree of nondiabaticity of 
the ionization. The peak field strength is Fo = 0.2 a.u., and 
the pulse duration is 400 a.u. The corresponding Keldysh 
parameter 7 is shown for different regions. 



low energies the total ionization probability is reproduced 
exactly by considering only the diabatic state (region I). 
For higher energies around 1 eV (region II), the differ- 
ence increases significantly, indicating that nondiabatic 
effects start to become important. 



2jt 



d ¥> r[/(i) costp], 



(19) 



where T[F] is the instantaneous ionization rate. Hence, 
the fact that the ADK theory of tunneling ionization 
and similar approaches cannot reproduce the correct (di- 
abatic) ionization rate for fe w-cy cle pulses is not due to 
coupling to higher states [24|, (25| , but rather because the 
pulse envelope changes dramatically within one cycle. In 
this limit the rate cannot be averaged over one period as 
was done in Eq. (|19p . whereas for multi-cycle pulses it 
can be used in combination with Eq. (|18|) yielding 



P{t) w cxp 



dt' T[f(t') 



(20) 



Analyzing region I in Fig. [6] further, we observe that the 
ionization probability is not constant as a function of 
photon energy. But the population loss in region I is 
well described by the ionization out of \^ ). According 
to our argument above, the apparent frequency depen- 
dence is rather a dependence on the form of the pulse or 
analogously on the relation between the cycles and the 
pulse envelope, which appears in a pronounced way for 
few-cycle pulses. Preferably, to avoid confusion, it could 
be called form dependence. As we have seen, the ioniza- 
tion behavior for few-cycle pulses can be well understood 
from the dynamics of a single diabatic state. 



C. Nondiabaticity and the special case of few-cycle 
pulses 

In order to find a common way of speaking we incorpo- 
rate the Keldysh parameter in our considerations, which 
has been used as an adiabaticity parameter. Following 
our language of the adiabatic representation, the ioniza- 
tion in the tunneling regime, 7 <C 1, is diabatic rather 
than adiabatic. We conclude from Fig. [6] that in the 
region where 7 « 1 the relative difference between the 
results calculated from the diabatic ionization rate via 
Eq. (|18p and from the solution of the Schrodinger equa- 
tion is greater than 10%. This is a clear sign of non- 
diabatic behavior. Already for 7 « 0.17 the diabatic 
ionization probability starts to differ slightly from the 
total ionization probability. For a fixed pulse duration 
we can also divide the frequency range according to the 
number of cycles in the pulse. Starting from the highest 
frequencies studied here we have multi-cycle pulses, until 
we reach few-cycle pulses at a photon energy of « 0.8 eV. 

The dynamics for few-cycle pulses is commonly con- 
sidered to be nonadiabatic (in our language this trans- 
lates to nondiabatic) [12J, |25j . We find that even for few- 
cycle pulses the tunneling is completely diabatic. In the 
framework of ADK theory and other approaches [5l[ the 
ionization rate T(t) is obtained by integrating over one 



V. CONCLUSION 

We have studied the dynamics of tunneling ionization 
in atoms and have found that, within the framework of 
the adiabatic representation, it is diabatic rather than 
adiabatic. We have identified two distinct ionization 
regimes depending on their diabatic behavior. In particu- 
lar we have characterized the transition from the diabatic 
to the nondiabatic regime. 

In the low-frequency limit the total ionization proba- 
bility is reproduced by the contribution of the tunneling 
probability of one single diabatic state. This means that 
in this regime there are no significant transitions to other 
diabatic states. For few-cycle pulses, the ionization prob- 
ability depends on the frequency for a fixed pulse dura- 
tion. However, this is not a nondiabatic effect, but the 
effect stems from the form dependence of the pulse, and 
the consequent fact that the rate cannot be averaged any 
longer over one period. 

When nondiabatic transitions start to happen, the dif- 
ference between the diabatic state ionization probability 
and the total probability increases dramatically. For fre- 
quencies in the range of the binding energy of the atom 
one-photon absorption can occur which is a completely 
nonadiabatic and even nondiabatic process. Already for 
parameters 7 ss 0.17 the diabatic ionization probability 
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starts to differ noticeably from the total ionization prob- 
ability, even though the pcrturbative multiphoton regime 
is not yet entered. From the perspective of the adiabatic 
representation, the Keldysh parameter is found to be an 
approximate measure of diabaticity. 
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Appendix A: Propagation methods 

This short section provides supplementary details on 
the technical features involved in our calculations. We 
need to find the exact solution of the Schrodingcr equa- 
tion [sec Eq. (JTJ] for the time-dependent Hamiltonian 
[see Eq. ©]. Via the propagation of the wavefunction 
in time, where the full dynamics are considered (adia- 
batic and nonadiabatic), we obtain the numerically ex- 
act solution. In the adiabatic approach the instantaneous 



Hamiltonian for different field strengths (corresponding 
to different time steps) is diagonalized. In the latter case 
the adiabatic eigenstates are obtained, so that the adia- 
batic dynamics are studied. 

For both methods we employ the time-dependent con- 
figuration interaction singles (TDCIS) schem e 1 53 . [53| . 
Starting with the Hartree-Fock ground state [54[ as the 
field- free ground state, |$o)i we include one-particle-one- 
hole configurations |55| , describing the excitations of the 
system. The TDCIS wavefunction reads: 



|*(t)) = a (t)|$o(*)> 






aUWm), (Al) 



where the index i symbolizes an initially occupied or- 
bital, whereas a denotes a virtual orbital to which the 
particle can be excited. This means that we consider 
only configurations where one particle is singly excited, 
thereby creating one hole. We use the software packages 



ARPACK [56j and LAPACK [57| to calculate the adi- 
abatic eigenstates. The two methods used permit us to 
compare the explicitly adiabatic dynamics with the exact 
calculation, where all dynamics arc included. 
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